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We carry out a covariant calculation of the measurable relativistic effects in an orbiting gyroscope 
experiment. The experiment, currently known as Gravity Probe B, compares the spin directions of 
an array of spinning gyroscopes with the optical axis of a telescope, all housed in a spacecraft that 
rolls about the optical axis. The spacecraft is steered so that the telescope always points toward a 
known guide star. We calculate the variation in the spin directions relative to readout loops rigidly 
fixed in the spacecraft, and express the variations in terms of quantities that can be measured, to 
sufficient accuracy, using an Earth-centered coordinate system. The measurable effects include the 
aberration of starlight, the geodetic precession caused by space curvature, the frame-dragging effect 
caused by the rotation of the Earth and the deflection of light by the Sun. 



I. INTRODUCTION 



Gravity Probe B - the "gyroscope experiment" - is a 
NASA space experiment designed to measure the gen- 
eral relativistic effect known as the dragging of inertial 
frames. The experiment will place into Earth orbit a 
spacecraft containing four gyroscopes and a telescope (for 
details of the experiment see fjL 0|) • The gyroscopes are 
4 cm diameter fused quartz or single crystal silicon balls, 
machined to be as spherical as possible (to tolerances 
better than one ppm), and coated with a thin film of 
niobium. At the low temperatures provided by a dewar 
of superfluid Helium surrounding the gyroscope assem- 
bly, the niobium is superconducting, and each spinning 
gyroscope has a London magnetic moment parallel to its 
spin axis (stray magnetic fields, such as those due to the 
Earth, will have been suppressed by many orders of mag- 
nitude before the niobium goes superconducting). The 
orientations of the gyros' magnetic moments are mea- 
sured via changes in the magnetic flux through supercon- 
ducting current loops that encircle each gyro, and that 
are attached to the gyroscope housing. The four gyros 
will initially have their spins aligned parallel to the sym- 
metry axis of the spacecraft, which coincides with the 
optical axis of the telescope mounted on the end of the 
spacecraft. In order to average out numerous unwanted 
torques on the gyros, the spacecraft will rotate about its 
symmetry axis at a rate of between one and 10 times per 
minute. 

According to general relativity, a perfect gyroscope in 
orbit about the Earth will precess relative to distant in- 
ertial frames because of two effects. The first, and most 
important, is the dragging of inertial frames, caused by 
the rotation of the Earth, also called the Lense-Thirring 
effect. Over a time of one year, for a polar orbit of about 
640 km altitude, this causes the gyroscopes to precess 
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in an East- West direction by around 42 milliarcseconds 
(mas). The second is the geodetic precession, caused by 
the curved spacetime around the Earth. This effect re- 
sults in an annual precession in the North-South direction 
of about 6600 milliarcseconds. 

The reference direction to the "distant stars" is pro- 
vided by the on-board telescope, which is to be trained 
on a star IM Pegasus (HR 8703) in our galaxy. One im- 
portant feature of this star is that it is also a strong radio 
source, so that its direction and proper motion relative to 
the larger system of astronomical reference frames can be 
measured accurately using Very Long Baseline Radio In- 
terferometry (VLBI) . In order to reduce telescope errors 
and biases, the spacecraft will be "steered" by attitude 
control forces to keep the centroid of the stellar image 
centered on the telescope axis. 

However, because the spacecraft orbits the Earth and 
the Earth orbits the Sun, the apparent direction to the 
star will vary in a periodic manner because of stellar aber- 
ration, which amounts to five arcseconds from the orbital 
motion, and to 20 arcseconds from the annual motion of 
the Earth. These are both substantially larger than the 
frame-dragging effect being measured. But, far from be- 
ing an annoying ancillary effect getting in the way of 
testing relativity, the aberration of starlight is in fact 
central to the success of the experiment ||. The reason 
is that the SQUID magnetometers that are connected to 
the current loops surrounding the gyros do not measure 
angle, they measure voltage, induced by the varying mag- 
netic flux threading the loops. They must therefore be 
calibrated to provide a scale factor that converts voltage 
to angle. Because the orbit of the spacecraft, the loca- 
tion of the star and the speed of light can all be estab- 
lished to high precision, the aberration of starlight can 
be predicted accurately. It can also be separated from 
the relativistic precessions because of its unique tempo- 
ral variability over the 12+ months of the science phase 
of the mission. Through this means, the experiment can 
be calibrated with an accuracy sufficient to achieve the 
desired overall accuracy. 

In such a complex and delicate experiment, there are 
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numerous sources of error, including anomalous, non- 
relativistic torques on the gyros, variations in electronics 
responses, errors in guide star position and proper mo- 
tion, effects of accumulated charges on the rotors from 
cosmic rays and other charged particles, and so on. Nev- 
ertheless, the team that is mounting this experiment 
- Stanford University, Lockheed-Martin Space Systems, 
and NASA - has proposed that a measurement of the rel- 
ativistic precessions at the level of 0.4 mas/yr is feasible; 
this would provide a test of relativistic frame dragging 
at the one percent level and of the geodetic effect at the 
level of 6 parts in 10 5 . 

The precession of gyroscopes according to general rela- 
tivity was first calculated independently by Pugh and by 
Schiff |Q, H ^| , who calculated the precession relative to 
a given inertial coordinate system. Wilkins |tJ analysed 
the precession in more detail, attempting in particular 
to ascertain how much of the geodetic effect could be re- 
garded as Thomas precession and how much was due to 
the curvature of space. He was also the first to attempt 
a more-or-less coordinate-free calculation by comparing 
the spin directions with the directions of incoming light 
rays from distant sources. 

Recent discussions of gyroscope precession generally 
fall into three groups. One group provides formal calcu- 
lations of the precession of a gyroscope relative to a local 
tetrad or congruence of geodesies, sometimes with ap- 
plications to specific spacetimcs such as Schwarzschild, 
Kerr, or de Sitter §, g, ||, 0, |l|, but seldom refers 
to specific aspects of Gravity Probe-B. Another group 
makes coordinate-dependent calculations of specific pre- 
cession terms (generally in the post-Newtonian approxi- 
mation) , and discusses details of the various terms, fre- 
quently with reference to Gravity Probe-B jl3| |lj, [D| 
[L6|. This group also includes unpublished calculations 
by the GP-B team in the course of developing the data 
analysis procedures for the experiment. A third group 
mixes calculations of precessions referred to a formal 
tetrad (sometimes including comparisons with the di- 
rection to a distant star) with applications to Gravity 
Probe-B; these include the standard textbook presenta- 
tions ]l7], |l8|, |2^, plfl , as well as papers designed in 
part to sort out the meaning of the various precession 
terms - geodetic, frame-dragging and Thomas precession 

ill- 

To date, we know of no calculation that, on the one 
hand, takes a fully covariant approach, and that on the 
other hand, incorporates critical operational aspects of 
the GP-B experiment, and performs an end-to-end cal- 
culation in terms of observable quantities. Such a co- 
variant calculation is of interest for several reasons. It 
would verify that this specific spacecraft experiment is 
actually measuring the observable general relativistic ef- 
fects it claims to be measuring. In addition, the effects of 
aberration are of order w/c and are not relativistic (they 
depend only on the finite propagation velocity of light), 
however they are subject to relativistic corrections which 
are of order (v/c) 2 , comparable to the two general rel- 



ativistic effects being measured. They should therefore 
be accounted for in a manner that does not depend on 
coordinates. 

We have carried out such a calculation. In Section II, 
we set up an orthonormal basis of vectors tied to the 
rolling spacecraft and use them to characterize the mea- 
surable components of the gyroscope spins as read by 
current loops and SQUIDs on board. Section III deals 
with the incoming light from the guide star and incorpo- 
rates the effects of spacecraft pointing on the spacecraft 
basis vectors. In Section IV we calculate the aberration 
and relativistic effects on the observed spin components 
using the parametrized post-Newtonian framework in an 
Earth- fixed reference frame. Section IV discusses the var- 
ious effects and estimates their magnitudes, and Section 
V g ives concluding remarks. We use the conventions of 
[tl9| and units in which G = c = 1. 



II. SPACECRAFT BASIS VECTORS AND 
GYROSCOPE SPINS 

We consider a spacecraft with four- velocity u, whose 
orientation is described by three spacelike basis vectors 
orthogonal to u that are rigidly tied to the spacecraft 
p3[ . The vector £ points along the spacecraft symme- 
try axis, while p and rj are orthogonal to £ and to each 
other. These four vectors e Q form an orthonormal basis, 
with e a ■ ep = r] a p, where rj a p is the Minkowski metric, 
valid in the local inertial frame of the orbiting space- 
craft. The spacecraft moves on a geodesic (we ignore at- 
mospheric drag and translational attitude control forces), 
with Vuu = 0. Because the spacecraft is rotating about 
its symmetry axis, it acts as a gyroscope, so that its spin 
axis £ is parallel transported, except for attitude-control 
forces that keep the telescope pointed toward the guide 
star, so that Vu£ = t, where t denotes the effect of those 
forces on £ ^4|. This, together with the orthonormality 
of the basis and the geodesic equation, are sufficient to 
establish that 

V U T7 = up - (*7-t)£, 

V u p = -uiri-(p-i)t, (I) 

where 

U) = p V U J7 = V V u p (2) 

is the locally measured roll angular frequency of the 
spacecraft. Defining the spacecraft roll phase 9 = J ujdr, 
where r is spacecraft proper time, we define "fixed" roll 
reference axes which do not rotate with the spacecraft: 

rj = r\ cos 9 — p sin 9 , 

p = r\ sin6> + pcos 9 . (3) 

The initial phase is chosen so that one of the vectors 
(t7 ) points orthogonal to the plane of the polar orbit, or 
in the East- West (EW) direction, while the other (p ) lies 
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Guide star 




FIG. 1: Basis four-vectors for Gravity Probe-B 

in the orbital plane, or the North-South (NS) direction 
(see Fig. 1). The orbital plane is chosen so that the guide 
star, and hence the vector lie in the plane. (The guide 
star HR 8703 is at declination 16.84° and right ascension 
343.26°.) Then one can show that these roll-fixed vectors 
evolve according to 



VuPo = -(Po-t)£> 
V u £ = t, 



(4) 



where we have used the fact that t • u = 0. 

Each gyroscope is characterized by a spin four- vector S 
that is purely spatial in the spacecraft frame (S • u = 0); 
this amounts to a specific choice of its center of mass. 
We assume that gyro torques are sufficiently small that 
we can treat the gyros as perfect, so that S is parallel 
transported (VyS = 0). Then the rate of change of 
each component of S on the roll-fixed basis is given by 
d(S ■ e a )/dr = S • Vue tt ; combining with Eqs. (||) we 
obtain the equation of motion for the spin components, 



Sr\ — 
Spa = 
St = 



—trio S% > 
—t po , 

trjo Srjo + tpo Sp j 



(5) 



where S Vo = S • r/ S Po = S • p , and = S • £, with 
parallel definitions for i^, etc., and where overdot de- 
notes d/dr. In the absence of torques on the spacecraft 
(t = 0), S vo = S Pa = — 0, as expected, since the gyro 
spins and spacecraft basis vectors will undergo identical 
transports through spacetime. 



III. TELESCOPE POINTING 

The spacecraft is oriented so that the optical axis of the 
telescope is aligned parallel to the direction of an incom- 
ing photon from the guide star. Let A be the null tangent 
vector to the incoming photon satisfying the geodesic 
equation V^A = with A • A = 0. In the spacecraft 
frame, the direction of the incoming photon is given by 
a unit normalized spatial vector which is the projection 
of A onto the spacecraft basis, namely 



(6) 



(-U- A) 



Then, using Eqs. (Q) the equations of motion for the 
roll-fixed components of A are given by 



(d/dT)X Vo 
(d/dr)X Po 
(d/dr)X 6 



e 1/2 (VuA)-77 -^ A c , 
e 1/2 (V u A).p -f po A ? , 



where we have introduced the small parameter e 1 / 2 ~ 
(v/c) to characterize the leading effects on VyA, namely 
aberration. 

In the operation of the experiment, attitude control 
forces orient the spacecraft so that the telescope always 
points toward the apparent location of the guide star, i.e. 
so that A = — £, or so that A p = A-77 = 0. This will be 
only approximately true because of pointing errors and 
occasional deliberate "dithering" of the spacecraft for op- 
erational and other reasons; these errors are expected to 
be on the level of a few to tens of mas. We introduce the 
small parameter 8 ~ 10~ 7 to characterize these pointing 
errors. Then to drive A Po and X Vo and their time deriva- 
tives to zero, we require attitude control torques given, 
from Eqs. @, by 

t no = e^AT^VuA^o + ^At^), 



c 1 / 2 Ar 1 (V u A)-p +*(At P0 ) ! 



(8) 



where At m and At po are the roll-fixed components of the 
order-<5 residual "pointing error" and dithering torques. 
Substituting Eqs. (||) back into Eqs. (Q), and using 
the fact that A is a normalized vector, we find that 
(d/dT)\ m = -5(M no )\£, (d/dT)X Po = -6(At Po )X^, and 
(d/dr)X^ = S(At Vo X Vo + At Po X Po ). From these equations 
we obtain 

X V0 (t) = 0(6) , X po (r) = 0(5) , A 5 (r) = -l + 0(<5 2 ) . (9) 
Then substituting Eqs. (§|) into (|), we obtain 

S V0 = -e 1 / 2 (S i /X i )(V u \)- Vo -6(At m )S i , 
S P0 = -e 1 / 2 (S 6 /X € )(V u \)-Po-S(At Po )S 6 , 



eV2[S no (V u A)- 77 + ^o(V,jA)-p ]A, 



-l 



-S(At Vo S Vo 



At Po S Po ) 



(10) 
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But since S Vo ~ S Po ~ e 1 / 2 + (5, we have ~ e + e 1 / 2 ^ + 
5 2 + . . ., and so, choosing the vector S to be normalized to 
unity (under parallel transport it has constant norm), we 
can write = 1 + 0(e,e 1 ^ 2 S,5 2 ). Substituting into Eq. 
( |l0| ) , we obtain our final equation for the evolution of the 
gyro spin components relative to the roll-fixed spacecraft 



e 1/2 (V u A).r, -<5(At ))0 ) 

+0(e^ 2 ,eS,e 1 / 2 S 2 ), 

e 1/2 (V u A).p -<5(Ai Po ) 

+0(e^ 2 ,eS,e^ 2 S 2 ). 



(11) 



These physically measurable components S„ and S po 
translate directly into outputs of the SQUID magnetome- 
ters. The expressions are generally covariant, and thus 
may be evaluated in any chosen coordinate system. 

Interestingly, in this picture all the relativistic effects 
arise mathematically from the covariant derivative of the 
incoming photon world line along the spacecraft four- 
velocity, and not in the precessions of the gyroscope spins. 
This is because we have chosen the rolling spacecraft it- 
self to define our basis; this was a natural choice because 
all the measuring instruments (current loops, SQUIDs) 
are rigidly fixed to the spacecraft. But since a rolling 
spacecraft is itself a gyroscope, it and the superconduct- 
ing gyros precess together in the spacetime around the 
Earth. Hence, in this basis, the relativistic effects ap- 
pear to originate elsewhere. The superconducting gyros 
are still crucial, of course, because, by design, they are 
closer to being perfect gyros, and thus are less affected 
by anomalous torques. Here the spacecraft plays a role 
similar to that of an oscillator clock that is slaved to an 
ultraprecise atomic standard (such as a hydrogen maser, 
or atom fountain clock); the oscillator may drift or be 
deliberately steered from time to time, but its phase can 
always be linked precisely back to the underlying stan- 
dard. Similarly, while the spacecraft "gyro" is steered 
to follow the guide star as the star's apparent location 
varies because of aberration and relativistic effects, the 
spacecraft axis direction is precisely linked to the super- 
conducting gyros via the output of the SQUIDS. 

We now turn to the relativistic and aberration effects 
that are being measured. 



IV. POST-NEWTONIAN CALCULATION OF 
MEASURABLE RELATIVISTIC EFFECTS 

We evaluate the quantity VyA in the post-Newtonian 
approximation using the parametrized post-Newtonian 
(PPN) framework jljj]. The metric of spacetime in the 
solar-system rest-frame is taken to have the form 

5 oo = -l + 2U + 0(e 2 ), 

g 0l = -(2 + 2 7 +ia 1 )U l +0(e 5 / 2 ), 

9ij = (l + 2 1 U)S. l] +0(e 2 ), 



where 7 and a\ arc PPN parameters whose values depend 
on the theory of gravity, and where we have restricted at- 
tention to theories of gravity with suitable global conser- 
vation laws, and have ignored "preferred- frame" terms in 
the metric (see [jl9|, §9.1 for discussion of these effects). 
We have used a gauge analogous to harmonic gauge, in- 
stead of the conventional PPN gauge (see [jl9| §4.2 for 
discussion). For an orbiting, rotating, nearly spherical 
Earth, and a static, spherically symmetric Sun, the po- 
tentials U and V have the form 



U = Mq / r + M<$ I R , 

V = M ffl v e /fl-lxxJ ffi /i? 3 : 



(13) 



where Mq is the mass of the Sun, M®, v® and J 9 are 
the mass, velocity and angular momentum of the Earth; 
X = x— Xq and R = |X| denote the satellite-Earth vector 
and distance, while x and r = |x| denote the satellite-Sun 
vector and distance. In fact, the quadrupole moment of 
the Earth has a measurable effect on the geodetic pre- 
cession and must be included; as this has been treated 
thoroughly elsewhere 0, we will ignore it here. 

From the geodesic equation for photons, d\ l /dt + 

( T U ~ T % Xl ) X,iy = °> whcrc AM = d ^/dt, we obtain 
the solution (see [[L9) §7.1 for details): 



A 



1, 
n'[l 



(l + 7 )C/]+2> i + 0(e 3 / 2 ) 



(14) 



where n l is a Cartesian unit vector describing the initial 
direction of the photon far from the Earth in the global 
PPN coordinate system, and V 1 represents the deflection 
of the ray, obtained by solving the equation 



~~dt 



(l + 7 )[V 4 [/-n l (n- VZ7)]. (15) 



along n 1 



(12) 



The unit vector n points from the source toward the 
spacecraft, and hence may be written n = no + x • Vno + 
. . ., where no is a constant unit vector toward the solar- 
system barycenter and where \7 l n J = (i5 y — n})n Q )/D, 
where D is the distance from the solar-system barycen- 
ter to the guide star. The correction term accounts for 
the parallax of the guide star; for HR 8703 it amounts to 
about 10 mas ~3x 10~ 7 . Because it is so small, we will 
treat it as effectively a pointing error term, of order S. 

Letting the spacecraft have ordinary velocity v l — 
u l /u°, where u° = 1 + v 2 /2 + U + 0(e 2 ), we find that the 
components of A are given by 

A = n • v - [v 2 - (n • v) 2 ] + 0(e 3 / 2 ) , 
A* = n l - [1/ - n l (n ■ v)] + n l [(n ■ v) 2 - v 2 /2 - jU] 
+V l + 0{e 3/2 ) . (16) 

To calculate the covariant quantities (VyA) • r) and 
(VuA) • p in Eqs. (|TT|) , it suffices to evaluate them in 
the PPN coordinate system. We first evaluate (VuA) a = 
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it^Afg. Calculating the Christoffel symbols from the 
metric (|l^) to the needed order, and using the fact 
that dU/dt = dU/dt + v • VU ", and, by virtue of the 
equation of motion dv/dt — VU, that (v • n)Vf/ = 
l/2{d[v(v • n)]/dt + n x (VZ7 x v)}, we find 



(V U A)° 
(V U A) 4 



0(vS/U) 
d 



dt 



(x-VK-v l +n J (n-v) 



+n 4 {(n-v) 2 -i« 2 }-^(n-v) 

+ i(2 7 +l)[nx (v x VU)f 

+ (l + 7+Jai)[nx (Vx V)] 4 
+0(UVU) , 



(17) 



where, to the order needed, we can convert from space- 
craft proper time r to PPN coordinate time t. 

We also need to determine the time dependence of r) n 
and p Q ; even though they are roll-fixed vectors, they 
will vary slightly because of pointing of the spacecraft 
to maintain alignment with the guide star. However, be- 
cause V U A is already of 0(e 1 ^ 2 ), we only need to evaluate 
r) and p to 0(e 1 ^ 2 ), which means accounting only for 
the effects of aberration-induced torques. Substituting 
Eqs. (g) into (Q), and using the fact that £, p , r] , and 
u form an orthonormal tetrad, we obtain 



V u r? 



e 1/2 V U A • Vo t + 0(6) , 
e 1/2 V U A • p | + 0(5) , 
e 1/2 ( V U A + V U A uu 
+0(8) . 



(18) 



We use the fact that 

\/3 _ ^ 



tors, (V u e (ct) )' 3 = d& [a) /dT 
Christoffel symbols and Eqs. 
tain 



for each of these three vec- 
' substitute the 



1 - 5 e (a) U 



), and integrate, to ob- 



( m f(t) = v.fja+0(e,S), 

M l (t) = (fj y- v .fj C+O(e,6), 

( Po )°(t) = vp + O(e,<5), 

(po)*(t) = (p Q f-v-p Q e + 0(e,5), 

e(t) = v£ + 0(e,5), 

C(t) = ? + v *-? y .Z + 0(e,8), 



(19) 



where fjo, po and £ are spatial 3- vectors in the PPN frame 
that describe the orientation of the spacecraft at a cho- 
sen initial moment of time. They are orthonormal in the 
Cartesian sense to 0(e). We have used the fact, obtain- 
able from Eq. (|), that n = -f + 0(e, 5), n-fj = 0(e, 5) 
and n • p = 0(e, 8). 

Substituting Eqs. @ and @ into (|ll]), converting 
from d/dr to d/dt using the property that t = £[l + 0(e)], 
and integrating with respect to t, we obtain, finally 



So 



+S(D- 1 x -fj -% 10 ) 
+0(e^ 2 ,eS,e i / 2 S 2 ), 
SM-e^vpo 

+e Q v ' Pov ■ f + U ■ p + V ■ po 

+^- 1 x-p -r po ) 

+0(e 3 / 2 ,e8,e^ 2 5 2 ). 



(20) 



The first term in Eqs. ( |20| ) is the initial misalignment of 
the spin with the spacecraft axis, which is expected to 
be smaller than one arcsecond. The next term, of order 
e 1 / 2 , and the first of the O(e) terms, are the aberration 
term and its relativistic correction. The term involving 
TL is the integrated relativistic precession projected onto 
?7o and po, given by Tt = J 1 (dH/dt)dt, where dJijdt is 
the precession rate given by 

^ = -i(2 7 +lKx(vxVC/)-(l+ 7 + ^ 1 )fx(VxV) . 

(21) 

The terms T> ■ t]q and V ■ po in Eqs. (^0|) denote the con- 
tribution of the deflection of light. The 0(6) terms T Vo 
and T po denote the integrated angular offset caused by 
pointing errors, while the terms D _1 x • fjo and _D _1 x • po 
denote the effects of parallax. The spatial vectors defined 
in Eqs. (|2(]) are referred to the solar-system- fixed PPN 
coordinate system. However, to the order needed, we 
can write v = v ffi + v s + 0(e 3 / 2 ), where v ffi is the Earth's 
orbital velocity relative to the Sun, and v s is the space- 
craft's velocity relative to the Earth. Furthermore, the 
components of the various 3-vectors in the solar-system 
basis are equal to the respective components in a non- 
rotating Earth-comoving basis, apart from corrections of 
order e; hence to the order needed in Eqs. (|20|) we can 
treat all spatial vectors as defined in the geocentric basis. 



DISCUSSION OF RELATIVISTIC EFFECTS 



s vo (o) 



'Vn 



The first non-constant term in Eqs. ( |20D is the aberra- 
tion of starlight, consisting of both an annual aberration 
with amplitude of v ffi ~ 20 arcseconds, and an orbital 
aberration with amplitude v s ~ 5 arcseconds. But be- 
cause the position of the guide star can be determined 
from VLBI and the orbit of the earth and spacecraft can 
be determined from the ephemeris and GPS to high ac- 
curacy, these terms may be predicted to an accuracy bet- 
ter than that needed to measure the relativistic effects. 
Furthermore, their variation with time over a year can 
be separated from the variation of the relativistic terms, 
hence the aberration can be used to determine the "scale 
factor" that relates the output of the SQUIDs (the ac- 
tual data) to the desired angles S Vo and S Po . The sec- 
ond term in Eqs. (p0|) is the relativistic correction to 
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aberration, whose orbital part dominates, with ampli- 
tude sin/3 ~ 0.4 mas (/3 is the ecliptic latitude of the 
guide star), which is comparable to the accuracy goal of 
the experiment; this term must be included in the data 
analysis pfjj. Terms of order Vq,v s and v 2 and terms cubic 
in v can be ignored. 

For the relativistic precession terms, dH/dt can be sep- 
arated into various contributions using Eqs. (|l3|): 



dH 
~db 



i? 3 

._ ( 2 7 + 1) _o Va 



i(27+l)|fv s xX 



■5(7+1 



1 s 



J© - 3NJ© ■ N 

i? 3 



x£, (22) 



where N = X/i?. Since 7i and £ both refer to the di- 
rection of the spins, this equation has the general form 
dS/dt = CI x S, so that the quantity in braces in Eq. (|2^ ) 
can be regarded as the precession rate vector. 

The first two terms produce periodic contributions to 
the signal at the frequencies u) s ±ui($, where w s and Ct>© <C 
u) s axe the orbital angular frequencies of the spacecraft 
and Earth, respectively. Integrated over time, these con- 
tribute periodic terms at approximately the orbital fre- 
quency with amplitudes of v^M^/2R 2 uj s ~ v^v s /2 ~ 0.3 
mas and w s M Q /r 2 a; s ~ v^(R/r) ~ 10~ 4 mas, respec- 
tively. Averaged over time, these are negligible. 

The third term may be evaluated to sufficient accuracy 
at the center of the Earth rather than at the satellite; 
treating the Earth's orbit as circular yields a signal that 
grows at the constant rate of (3/2)u^W0(h x £~) ~ 19 
mas/yr in a direction perpendicular to h, the normal to 
the ecliptic plane 0, ^?J. This term is the solar geodetic 
effect (sometimes called the dc Sitter precession); it is 
responsible for an analogous precession in the lunar orbit 
that has been measured to around 0.7 % using Lunar 
laser ranging f28|| . 

The fourth term is the main geodetic effect, at the 
constant rate of (3/2)VgW s (h' x £) ~ 6600 mas/yr, in a 
direction perpendicular to h', the normal to the orbital 
plane. Note that if 770 is chosen to be the EW basis vector, 
which is perpendicular to the polar orbit, then the main 
geodetic effect does not appear in the output S Vo , in other 
words it produces a purely S P(n or NS signal. There is 
an additional contribution to this precession arising from 
the Earth's oblateness at about 7 mas/yr, which must be 
taken into account jj], [l5|, [la |25|, |29| , |3(J 

The final term in Eq. (|22j)is the Lense-Thirring or 
frame-dragging effect. It has both a constant term and a 
periodic term at twice the orbital frequency, but for a cir- 
cular polar orbit, it produces an orbit-averaged constant 
rate of J© x t;/2R 3 = J©cos(5/2i? 3 ~ 41cos<5 mas/yr, 
where S is the declination of the guide star. If the plane 
of the polar orbit is oriented so that the guide star loca- 
tion is in the orbital plane, then the signal is purely S m , 
or EW. 

The next term in Eqs. pel) is the deflection of light. 



For the ecliptic latitude of the guide star, this gives a 
maximum deflection caused by the Sun of around 20 mas, 
once in the one-year mission. This effect can be included 
as a term to be estimated in the data analysis - studies 
by the GP-B team indicate that it could be measured 
to about one percent - or it can calculated a priori us- 
ing general relativity and the known parameters of the 
Earth's orbit and the guide star location. Because the 
guide star will be occulted by the Earth, the Earth will 
also contribute a once-per-orbit deflection of around 0.5 
mas (which will compete with deflection caused by prop- 
agation of the signal through the Earth's atmosphere). 
Because it averages to zero, this effect can be ignored. 

These relativistic terms, geodetic precession, frame 
dragging and deflection of light, have been calculated us- 
ing the specific assumptions and constraints inherent in 
the PPN framework. An alternative, more phenomeno- 
logical viewpoint, would be to regard each as an indepen- 
dent effect with a separate phenomenological parameter 
to be measured in the data analysis. 

The remaining terms in Eqs. ( p0| ) are of order 8. The 
parallax term caused by the Earth's motion (around 10 
mas) must be taken into account, while that due to the 
satellite's orbit is negligible. Finally the pointing error 
terms will be controlled to better than 20 mas (rms), and 
in any event will be measured directly by the telescope; 
all that is needed then is a sufficiently accurate calibra- 
tion (via deliberate dithering of the spacecraft axis by 
known amounts) to convert telescope readout signals to 
the appropriate contributions T m and T po to S m and S Po . 



VI. CONCLUSIONS 

We have carried out a covariant calculation of the 
measurable output of an orbiting gyroscope experiment, 
and shown that the result may be expressed in terms 
of contributions that can be calculated consistently us- 
ing an Earth-based coordinate frame. The dominant 
contributions that are detectable within the stated 0.4 
mas/yr accuracy of the GB-P mission are: the aberra- 
tion of starlight and its special relativistic correction, the 
general relativistic geodetic precession (due both to the 
Earth and to the Sun) and the general relativistic frame- 
dragging. The results are in agreement with other, non- 
covariant calculations, including numerous unpublished 
calculations and simulations carried out by the GP-B 
team in the course of developing the data analysis sytem 
for the mission. Because all measurements are referred to 
instruments fixed on a rolling spacecraft, the measurable 
relativistic effects enter mathematically via their effect on 
the apparent direction of the light from the guide star. 
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APPENDIX: EFFECT OF ABERRATION ON 
MEASUREMENT OF ROLL PHASE 

The actual signals measured on the GP-B spacecraft 
are S v and S p ; these are then converted to roll-fixed sig- 
nals in the data analysis by forming the combination 

S Vo = S n cos 9 — S p sin 9 , 

S Po = S v sin 9 + S p cos 9 , (A.l) 

where 9 is the spacecraft roll phase. In practice this is de- 
termined using a pair of star-trackers fixed to the space- 
craft platform directed roughly perpendicular to the roll 
axis, that continually compare the spacecraft orientation 
to a library of known stellar positions, and provide a 
phase accurate to a few arcseconds. The angle 9 is the ac- 



cumulated phase in the onboard rotating frame since the 
initial state, so that any error in measuring 9 relative to 
its "true" value will generate corrections to the roll-fixed 
outputs given by AS m « —S Po A9, and A5 Po « S Vo A9. 

Consider a star-tracker whose optical axis is parallel 
to the rj vector. It will record a star when the incoming 
projected world line of the photon is anti-parallel to r\. 
i.e. when r\' = —A = n + v — nn ■ v. In the absence 
of spacecraft motion, we would have 77 = — n. If 9' is 
the measured roll phase and 9 is the "true" roll-phase 
(without aberration), then 

cos 9' — rj' ■ rj 

= cos 9 + v ■ psin6> , (A.2) 

hence A9 = v • p = v • po cos 9 — v • r]o sin 9. But since 
S Vo = S Vo (0) — e x / 2 v • 770 + 0(e, 6), with a parallel expres- 
sion for S Po , the product S^A^ is of magnitude at most 
10~ 8 ~ 2 mas, but most importantly is periodic at the 
roll frequency, and averages to zero. Aberration will be 
taken into account in determination of roll phase from 
the star trackers. 
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